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ABSTRACT 

We use the “Evolution and Assembly of GaLaxies and their Environments” (EAGLE) suite 
of hydrodynamical cosmological simulations to measure offsets between the centres of stel¬ 
lar and dark matter components of galaxies. We hnd that the vast majority (> 95%) of the 
simulated galaxies display an offset smaller than the gravitational softening length of the sim¬ 
ulations (Plummer-equivalent e = 700 pc), both for field galaxies and satellites in clusters and 
groups. We also hnd no systematic trailing or leading of the dark matter along a galaxy’s di¬ 
rection of motion. The offsets are consistent with being randomly drawn from a Maxwellian 
distribution with a ^ 196 pc. Since astrophysical effects produce no feasible analogues for 
the 1.62^Q'4g kpc offset recently observed in Abell 3827, the observational result is in tension 
with the collisionless cold dark matter model assumed in our simulations. 
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1 INTRODUCTION 

Observations of one galaxy in cluster Abell3827 (redshift z~0.1; 
Carrasco et al. 2010) revealed a surprising 1.62lg'4g kpc (68% CL) 
offset between its dark matter and stars (Massey et al. 2015). 
Such offsets are not observed in isolated field galaxies (Koop- 
mans et al. 2006; Gavazzi et al. 2007)^. However, offsets inside 
clusters are consistent with theoretical predictions from models of 
self-interacting dark matter (SIDM Spergel & Steinhardt 2000). As 
galaxies move through a cluster core, interactions with the cluster’s 
dark matter would create a friction and cause a galaxy’s dark mat¬ 
ter to lag slightly behind its stars (Massey et al. 2011; Kahlhoefer 
et al. 2014; Harvey et al. 2014), just like ram pressure causes gas 
to lag a long way behind stars in the Bullet Cluster (Clowe et al. 
2004, 2006; Randall et al. 2008; Harvey et al. 2015). Simple sim¬ 
ulations tailored to Abell 3827 support this prediction (Kahlhoefer 
et al. 2015, although current results operate under the limited as¬ 
sumption that the galaxy is on first infall). Many particle physics 
models of dark matter naturally predict low level self-interactions 
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^ A small number of galaxy quad lenses are not well-fitted by standard 
parametric models of dark matter centred on the optical emission. To fit lens 
RXS J1131, Claeskens et al. (2006) need to include a 0.044" offset or m = 
4 octupole term. With lens COSMOS J09593, Jackson (2008) achieved an 
acceptable goodness of fit only with a 0.063" offset and (an unrealistically 
large) external shear I7I = 0.25. However, in these isolated lenses the cause 
of these poor fits is more likely to be local substructure (Hezaveh et al. 
2013). An offset between mass and light would produce a relatively shallow 
core profile and possibly more detectable central images. Note also that the 
location of mass peaks is determined much more precisely by strong lensing 
than by weak lensing (George et al. 2012; Dietrich et al. 2012). 


(e.g. Tulin et al. 2013; Foot 2014; Boddy et al. 2014; Hochberg 
et al. 2014; Cline et al. 2014; Khoze & Ro 2014). If the interac¬ 
tion cross-section is considerable > O.lcm^/g, it could also re¬ 
solve small-scale issues in the predictions of inert, cold dark matter 
(CDM) models (see review by Weinberg et al. 2015). 

Cluster Abell 3827 was originally studied by Williams & Saha 
(201 1) because its light distribution is interesting, with the intention 
of developing a lens analysis algorithm but not with the expectation 
of measuring an offset (L. Williams 2015, pers. comm.). This is 
the only galaxy for which an offset has been detected, but it may 
also be the only galaxy in a cluster for which such a small offset 
could have been detected. The measurement requires three chance 
circumstances, each individually rare. 

• The cluster must gravitationally lens a well-aligned background 
galaxy with a complex morphology. The distribution of foreground 
dark matter (plus baryons) can be reconstructed from perturbations 
to this lensed image. 

• The cluster must contain a bright galaxy near the Einstein ra¬ 
dius. To enable precise measurements, it must intersect the lensed 
arcs and its mass must be a detectable fraction of the cluster. Since 
a single cD galaxy generally lies inside any Einstein rings, in prac¬ 
tice, this means a cluster with multiple cDs. 

• The cluster must be nearby, so small physical separations can be 
resolved. This reduces its efficiency as a gravitational lens. 

The interpretation of the observed offset in such radical terms 
as SIDM is clouded by the possibility of alternative explanations. 
First, gravitational lensing is sensitive to the total mass distribu¬ 
tion projected along the line of sight. The chance alignment of 
unrelated foreground/background structures has created apparently 
spurious features in other lens systems (Gray et al. 2001; Hoekstra 
2003; Host 2012). Similarly, source-lens degeneracies could lead to 
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a similar effect (e.g. Schneider & Sluse 2013, 2014). In Abell 3827, 
projection effects do not appear to be an issue: of the four galaxies 
at the centre of the cluster three have a total mass appropriate for 
their stellar mass, while the fourth (galaxy Nl) has a low mass at the 
location of the stars, but a similarly appropriate total mass slightly 
offset. Had this been a chance projection, there would be mass at 
the location of Nl (because of its own, non-offset dark matter) plus 
a second mass peak (and probably a luminous source). These are 
not seen. 

Second, a physical offset might arise even with collisionless 
dark matter, via the complex astrophysical processes operating in 
cluster core environments. Gas stripped from and trailing behind an 
infalling galaxy may self-gravitate and form new stars. This is not 
consistent with observations of Abell 3827, which has effectively 
zero star formation rate (Massey et al. 2015, Table 1). The differ¬ 
ent physical extent of dark matter and stars also leads to different 
dynamical friction, tidal gravitational forces, and relaxation times 
during mergers. Inside the complex distribution of Abell 3827, even 
normally linear effects like tidal forces could create or exacerbate 
small initial offsets. It could also be considered that the galaxy in 
question is undergoing one of two types of merger: 

• Coincidentally with the galaxy’s arrival near the cluster core , it 
has recently merged with a former satellite. The tightly-bound stars 
from the centre of the satellite have not yet had time to mix with 
the galaxy’s stars, and remain as a second peak randomly located 
within the total system. Simulated analogues of this are not consis¬ 
tent with observations, because the observed galaxy is best-fit in all 
bands by a single Sersic profile (Massey et al. 2015, Table 1). 

• The galaxy is about to merge with a more massive halo (the 
three more central galaxies of similar mass). In simulations, the 
dark matter from all the systems rapidly mix together into a single 
smooth halo. This is not consistent with observations, which still 
show the infalling galaxy’s dark matter, distinct from and further 
away from the other galaxies’ dark matter^. 

As a control test to determine whether more complex astro- 
physical effects could build an offset between galaxies and colli¬ 
sionless dark matter, we measure the 3D separation between galax¬ 
ies’ luminous and dark matter in the “Evolution and Assembly 
of GaLaxies and their Environments” (EAGLE) suite of hydrody- 
namical cosmological simulations (Schaye et al. 2015; Crain et al. 
2015). These simulations have been calibrated to match the masses 
and sizes of galaxies in the local Universe. The main EAGLE sim¬ 
ulation also reproduces the observed low-redshift luminosity func¬ 
tions (Trayford et al. 2015) and produces an evolution of the galaxy 
mass function in broad agreement with observations (Furlong et al. 
2015). The simulated galaxies display rotation curves in agree¬ 
ment with observations whilst the stellar and dark matter profiles 
of BCGs match observational data (Schaller et al. 2015a,b). Simi¬ 
larly, the tidal stripping and ram pressure stripping of the satellites 
as well as the AGN activity in the BCGs lead to a realistic popula¬ 
tion of galaxies in clusters (in terms of colours or SFR), indicating 
that the processes that could move matter around are broadly re¬ 
produced by the model. The EAGLE simulations are hence, an ideal 
test-bed to predict the relative positions of galaxies’ various com¬ 
ponents in a statistically meaningful way. 


^ Allowing a distinct dark matter peak for Nl fits the observations with 
X^/dof = 49.3/23, Bayesian evidence logjQ(iJ) = —26.4, and lens-plane 
(rmsi) = 0.26" (Massey et al. 2015, Table 3). A model without dark matter 
(but still stellar mass) is strongly disfavoured, with x^dof = 86.1/26, 
logj^pjE) = —100.7 and (rms^) = 0.34" (R. Massey 2015, pers. comm.). 


2 METHOD 

In this section we describe briefly the cosmological simulations we 
analysed and the method used to infer the centre of luminous and 
dark matter in galaxies. 

2.1 The simulation suite 

In our study, we use the main EAGLE simulation (Ref-L100N1504) 
and to explore field galaxies, clusters and groups, and the higher 
resolution simulation (Recal-L025N0752) to understand the con¬ 
vergence of our results. These cosmological simulations use a state- 
of-the-art treatment of smoothed particle hydrodynamics and set 
of subgrid models. The full description of the model is given in 
Schaye et al. (2015) and the rationale for its parametrisation is pre¬ 
sented in Crain et al. (2015); we only summarise here the aspects 
relevant to our study. The simulations assume collisionless dark 
matter, evolving in a flat ACDM cosmology with parameters from 
Planck2013 (Ade et al. 2014). The low (high) resolution initial con¬ 
ditions are generated at 2 = 127 using second-order Lagrangian 
perturbation theory in a 100® Mpc® (25® Mpc®) volume with a 
dark matter particle mass of 9.7 x 10® Mq (1.2 x 10® Mq) and 
initial gas particle mass of 1.8 x 10® Mq (2.2 x 10® Mq). The par¬ 
ticles are then evolved in time using the Gadget Tree-SPH code 
(Springel 2005). The Plummer-equivalent gravitational softening is 
set to e = 700 pc (e = 350 pc at higher resolution). 

The subgrid model in the EAGLE simulations includes 
element-by-element radiative cooling (Wiersma et al. 2009a), star 
formation obeying the Kennicutt-Schmidt relation (Schaye & Dalla 
Vecchia 2008), enrichment of the ISM via stellar mass loss 
(Wiersma et al. 2009b), feedback from star formation (Dalla Vec¬ 
chia & Schaye 2012), gas accretion onto super-massive black holes 
and the resulting AGN feedback (Booth & Schaye 2009; Rosas- 
Guevara et al. 2013). 

2.2 Identification of galaxies and their locations 

We find galaxies in the simulation via the SUBFIND algorithm 
(Springel et al. 2001). We identify all galaxies with stellar mass 
Mt > 10® Mq at 2 = 0, both in the field and in groups or clusters. 

We find the centre of galaxies’ matter distributions using an 
iterative ‘shrinking sphere’. We first identify all the star particles 
for each galaxy. We calculate their centre of mass and the distance 
of every particle to this centre. We then select only those parti¬ 
cles within 90% of the maximal distance to the centre of mass. 
Repeating this process, the search radius and the number of con¬ 
sidered particles decreases in subsequent iterations. This shrinking 
sphere procedure is repeated until the number of particles reaches 
200. This typically corresponds to a sphere of radius ~ 1 kpc, i.e. 
slightly larger than the softening length of the simulation. The cen¬ 
tre of mass of this final set of particles is considered to be the centre 
of the galaxy’s stellar distribution®. Similarly, we define the veloc¬ 
ity of the stellar distribution as the mass weighted velocity of the 
particles selected in the final iteration of the procedure. 

The same procedure is applied to each galaxy’s dark matter 

® As pointed out by Kahlhoefer et al. (2014, 2015), the choice of centroid- 
ing algorithm could produce varying results if dark matter does interact. Our 
identification of mass-weighted peaks in the stellar particles is both robust 
and the most comparable procedure to the identification of peaks in if-band 
luminosity-weighted observations (or other infrared bands in the absence of 
recent star formation). 
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Figure 1. The offset between the centre of the dark matter distribution and 
stellar distribution for galaxies with a stellar mass M* > 1O®M0. The dif¬ 
ferent panels correspond to field galaxies in the reference simulation (top 
left), the field galaxies in the higher resolution simulation (top right), satel¬ 
lites in clusters (bottom left) and in groups (bottom right). In each panel 
the vertical dashed line indicates the softening length used in the simula¬ 
tion. The arrows indicate the position of the 95% and 99.7% percentiles of 
each distribution. The offset seen is similar in field galaxies and clusters and 
is of the order of the softening scale of the simulation. Offsets larger than 
1.5 kpc correspond to fluctuations greater than 3cr. 

particles, to calculate the centre of their dark matter distribution. 
Finally, the offset between the dark and luminous component is 
defined as the distance between those two centres. We have verified 
that varying the minimum number of particles to define a galaxy 
centre from 100-500 and the shrinking ratio from 0.5-0.99 does not 
significantly affect our results. 


3 OFFSETS BETWEEN DARK MATTER AND STARS 

In the EAGLE Ref-L100N1504 simulation there are 12776 galax¬ 
ies with mass M* > 10® Mq, 1129 of which are satellites in 
clusters (halos with M 200 > 10^'* Mq), 3111 satellites in groups 
(halos with M 200 > 10^^ M©) and 7391 are field galaxies. 
The higher resolution Recal-L025N0752 simulation contains 618 
galaxies above our mass threshold. These four samples will be used 
to investigate environmental and resolution effects. 

3.1 3D offset between dark matter and stellar components 

The offsets between the centre of galaxies’ dark matter and their 
stars for our four sub-samples of galaxies is shown in Fig. 1. The 
distributions are consistent with being randomly sampled from a 
Maxwellian with distribution parameter o = 196 ± 2 pc (main sim¬ 
ulation) or cr = 126 ± 1 pc (high resolution simulation). Arrows 
indicate the position of the 95 and 99.7 (2 and 3cr) percentiles. In 
both cases, the typical scatter is smaller than the gravitational soft¬ 
ening length, indicated by a vertical dashed line. 


The distribution of offsets in the Ref-L100N1504 main simu¬ 
lation is remarkably similar for field galaxies (top left panel) and 
satellite galaxies in groups or clusters (bottom panels). This indi¬ 
cates that, at our resolution, the offsets are not influenced by envi¬ 
ronmental effects. Fewer than 5% of all galaxies display an offset 
larger than the gravitational softening length. Offsets larger than 

1 kpc are only found in 59 field galaxies (0.79%) and 17 satel¬ 
lites in groups and clusters (0.54%). Pushing these numbers to off¬ 
sets larger than 1.5 kpc, we find 15 field galaxies (0.20%) and 2 
satellites in groups and clusters (0.06%). A much larger sample of 
galaxies would, however, be required to characterise robustly the 
tail of the distribution. 

Offsets in the higher resolution Recal-L025N0752 simulation 
(top right panel) are smaller, with 95% of the galaxies displaying 
an offset smaller than 410 pc. Unfortunately, the smaller number 
of galaxies in that simulation volume does not allow for a thorough 
discussion of the position of larger percentiles. The results from this 
simulation indicate that the offsets seen in the main simulation are 
probably overestimated (at least for field galaxies) and that simula¬ 
tions run at a higher resolution (i.e. with a smaller softening length) 
would lead to galaxies with smaller offsets between dark matter 
and stars. However, the decrease in softening length by a factor of 

2 between our two simulations has only led to a decrease in me¬ 
dian offset by a factor of 1.5, indicating that even higher resolution 
simulations might not converge towards a negligible offset between 
components^. We nevertheless caution that the softening length is 
not the only scale setting the resolution of a simulation. Changes in 
the subgrid parameters and, sometimes, models between different 
simulations at different resolution are necessary to account for the 
newly resolved scales and have a non-trivial impact on the analysis 
of convergence. 


3.2 Offset along the direction of motion 

If the dark matter-stellar offset in Abell 3827 is due to SIDM, then 
not only will the centres of the galaxies and dark matter halos be 
offset but this offset should also be aligned with the direction of 
motion of the galaxy with the dark matter trailing the stars. Al¬ 
though the offsets observed in the EAGLE simulation are approach¬ 
ing the resolution limit set by the scale of gravitational softening, 
it is worth measuring whether the dark matter might be trailing or 
leading the galaxies in their motion. 

The offset between dark matter and stars, projected along the 
velocity vector of the stars, is shown in Fig. 2. In all four galaxy 
sub-samples, the distribution is symmetric and shows no bias to¬ 
wards leading or trailing motion of the dark matter. The distribu¬ 
tion and its mirror image are indistinguishable in a Kolmogorov- 
Smirnov (KS) test, with ap-value larger than 0.9. The scatters are 
(7 « 210 pc (main simulation) or cr « 128 pc (high resolution sim¬ 
ulation), in agreement with those for the Maxwellian 3D offsets. 
The dashed lines in the figure show Gaussians with the same mean 
and width as the measured distributions. The distributions of pro¬ 
jected offsets seem leptokurtic, with more offsets near 0 than the 


The offset of 300-400 pc found by Kuhlen et al. (2013) in high resolution 
zoom-in simulation of a single Milky Way-like galaxy is consistent with our 
findings. That offset from the centre of their dark matter distribution is ~3 
times larger than their softening, indicating that a small but non-zero offset 
might be found with sufficiently high resolution adopted in the simulations. 
We note, however, that the dark matter density profile of their galaxy is not 
monotonic; a result different from what is seen in other simulations. 
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Figure 2. The offset between the centre of the dark matter distribution and 
stellar distribution along the axis of motion of the stellar distribution for 
galaxies with a stellar mass M* > IO^Mq. The different panels corre¬ 
spond to the same subsets of galaxies as in Fig. 1. The dashed grey curves 
in the background show a Gaussian distribution with the same mean and 
standard deviation as the offset distributions. The distribution of offsets dis¬ 
plays no bias towards trailing or leading motion of the dark matter centre 
with respect to the luminous centre and only deviates from a normal distri¬ 
bution by displaying a positive kurtosis. 


equivalent Gaussian, but the small number of galaxies in the sim¬ 
ulation does not probe the tails of the distribution. The offsets are, 
thus, consistent with being randomly orientated and unaffected by 
the motion of the galaxy. Indeed, we find no preferred direction of 
offset, repeating the experiment by projecting the offset onto other 
axes like the dark matter velocity, direction to nearest neighbour, 
direction to cluster centre, etc. All offsets are consistent with ran¬ 
dom scatter. 


3.3 Detailed examination of satellite galaxies in the tail of the 
distribntion 

In our sample of satellite galaxies, we found 17 (2) objects out of 
3111 presenting an offset larger than 1 kpc (1.5 kpc). It is hence 
worth exploring whether these are just random fluctuations in the 
population or whether these larger offsets are seen as the result of 
an astrophysical process. The offsets of the 17 galaxies display no 
preferred direction with respect to the direction of motion, with a 
flat distribution of cos [6), where 6 is the angle of the offset from 
the velocity vector of the galaxy. 

The first of the two extreme outliers is a low mass (M* = 
3.9 X 10® Mq) extended galaxy (rso = 6.9 kpc). This galaxy is 
too diffuse at the resolution of the simulation for the centre-finding 
algorithm to return a sensible answer. Similarly, it would be diffi¬ 
cult to find the centre of the light distribution of a galaxy with such 
a flat profile in real observations. A galaxy like the one for which 
an offset has been observed in Abell 3827 is much more massive 


and less diffuse, making the presence of this specific outlier in our 
catalogue irrelevant for the scenario we are considering. 

The second extreme outlier is a giant elliptical galaxy with 
stellar mass M* = 1.5 x 10^^ Mq, located 130 kpc from the 
centre of its cluster. This galaxy experienced a recent merger with 
a smaller very concentrated satellite (M, ~ 2 x 10® Mq). 

The dark matter from the two galaxies has mixed, forming a 
smooth, virialised halo. The stars from the elliptical lie at the centre 
(within 200 pc) of this dark matter. However, the tightly-bound 
stars from the centre of the former satellite have not yet had time 
to mix with the stars from the elliptical. They instead remain as a 
peak in the outskirts of the stellar light distribution. This merger 
remnant is thus affecting the measurement of the peak of the light 
distribution but, at the time of measurement, it does not carry any 
dark matter. The large perceived offset is a temporary phenomenon 
due to the difference between the time taken to mix the stars in 
interacting galaxies and the time needed to mix their dark matter. 
This is the first merger scenario discussed in the Introduction, and 
is not consistent with details of the observations (except perhaps a 
short time window might exist during which distinct dark matter 
peaks still exist, but are offset from the light. This time window 
would make Abell 3827 even more rare.) 


3.4 Detailed examination of satellite galaxies in the cores of 
cinsters 

The simulation contains 50 (11) M* > x 10^® Mq satellite galax¬ 
ies within the central 100 kpc of groups (clusters). The statistics for 
this small sample are noisier, but they have a similar offset distri¬ 
bution and distribution of angles between offset and velocity vector 
as the full sample. The distribution of angles is consistent with uni¬ 
form and the distribution of offsets has a mean of 310 pc with a 
95 percentile at 690 pc, in remarkable agreement with the whole 
population. This sub-sample and the whole population are virtually 
indistinguishable in a KS test (p-value > 0.6). 

The closest non-BCG galaxy with M* > lO^^M© in the six 
simulated M 200 > 10^^ M© clusters are at clustercentric radii 
26 kpc, 92 kpc, 22 kpc, 58 kpc, 82 kpc and 54 kpc. These have off¬ 
sets between their stars and dark matter of 182 pc, 223 pc, 252 pc, 
198 pc, 320 pc and 284 pc, in apparently random directions. Look¬ 
ing in more detail at the two objects with the smallest clustercentric 
position, we find two elliptical galaxies of mass 1.5 x 10^® and 
4 X 1O^®M0 with low star formation rate and gas content. They 
both present an offset between their dark and luminous component 
smaller than 250 pc, unaligned with their direction of motion nor 
aligned with the radius to the centre of the cluster. We thus find no 
feasible analogues for Abell 3827 in the EAGLE simulation. 


4 SUMMARY & CONCLUSION 

Motivated by the measurement of a 1.62to'49 kpc offset between 
the stars and dark matter of a galaxy in Abell 3827 (Massey et al. 
2015), we investigated the relative location of these matter com¬ 
ponents in galaxies from the ACDM EAGLE simulation suite. Our 
results can be summarised as follows: 

(i) More than 95% of simulated galaxies have an offset between 
their stars and dark matter that is smaller than the simulation’s grav¬ 
itational softening length (e = 700 pc). The offsets are smaller still 
in our higher resolution simulation, indicating that our measured 
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values are likely upper limits. Even this state-of-the-art cosmolog¬ 
ical simulation has only just sufficient resolution to compare to the 
observations. 

(ii) Of the extreme objects with resolved offsets, fewer than 
0.54% (0.20%) of satellite galaxies in groups and clusters present 
a separation larger than 1 kpc (1.5 kpc). 

(iii) We find no systematic alignment between the direction of the 
offset and the direction of motion of the galaxies. Dark matter is 
statistically neither trailing nor leading the stars. 

(iv) We find no difference between field galaxies and satellite 
galaxies in groups and clusters. Astrophysical effects related to a 
galaxy’s local environment play no significant role in producing or 
enhancing offsets. 

(v) We find two types of outliers with extreme offsets: faint galax¬ 
ies for which the resolution of the simulations does not allow for the 
robust identification of a centre, and massive galaxies that have re¬ 
cently absorbed a smaller galaxy but haven’t yet mixed their stellar 
distributions. Neither of these outlier types match what is observed 
in Abell 3827. 

(vi) Looking specifically at the massive satellite galaxies close to 
cluster cores, we find no difference between these objects and the 
overall population of satellites or field galaxies. Environmental ef¬ 
fects do not seem to create offsets. 

Astrophysical effects, as modelled in the EAGLE simula¬ 
tion, produce no feasible analogue, for the galaxy observed in 
Abell 3827. Taking the best-fitting value for its observed offset, 
this galaxy would be a > 3 (t outlier in a ACDM universe with col¬ 
lisionless dark matter. Larger, higher-resolution simulations will, 
however, be needed to improve the sampling of the tail of the offset 
distribution and and to assess if the offsets measured in our simula¬ 
tion are biased high by limited numerical resolution. 

The observation is so far unique, and finding more systems 
in which similarly precise measurements can be obtained will be 
challenging. If more large offsets can be found and larger simula¬ 
tions confirm our findings, the case for an alternative dark matter 
model (e.g. SIDM) would be compelling. High resolution simula¬ 
tions including these models of dark matter would also be useful, to 
understand the processes that might have led to the observed offset 
in Abell 3827. 
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